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Abstract. We present a power spectrum analysis of interferometric 21 cm observations of the Ursa Major high- 
latitude cirrus, obtained with the Dominion Radio Astrophysical Observatory (DRAO) of Penticton (Canada). 
These high-resolution data reveal the intricate structure of the diffuse Galactic Hi, at angular scales from 1 
arcminute to 3 degrees. A filtering method based on a wavelet decomposition was used to enhance the signal-to- 
noise ratio of the data. The power spectra of the integrated emission and of the centroid velocity fields were used 
to deduce the three-dimensional (3D) spectral index of the density and velocity fields of the Hi cirrus. The spectral 
index is similar for the 3D density and velocity fields with a value of —3.6 ± 0.2. Using the Leiden/Dwingeloo 
observations, this analysis was extended to the whole North Celestial Loop (which includes the Ursa Major cirrus), 
showing that the scaling laws prevail from 0.1 to 25 pc. The centroid velocity and integrated emission fields show 
moderate correlation, with a maximum cross-correlation value of 0.44. 



1. Introduction 



It has been noted long ago ijChandrasekharl [l949; 
Ivon Weizsackerlll95l|l that the interstellar medium (ISM) 
has the dynamical properties of a turbulent flow. The high 
Reynolds number of the ISM has often been invoked to 
ascribe its self-similarity, observed with several tracers on 
scales from a few tenths to several hundreds of parsecs, to 
the presence of a fully developed turbulence. Nevertheless, 
despite the numerous efforts done to determine the dy- 
namical properties of interstellar gas, the properties of in- 
terstellar turbulence are still poorly known. This is mainly 
due to the difficulty to deduce three-dimensional prop- 
erties of the gas using projected quantities, but also to 
opacity effects and to the limited range of scales observed. 
Characterizing interstellar turbulence, by determining its 
statistical properties for instance, is still of major im- 
portance for the understanding of the structure and the 
kinematics of the ISM and its impact on the star forma- 
tion process. In particular the statistical properties of the 
ISM structure and kinematics can be used to identify re- 
gions with exceptional dynamical properties and study 
their impac t on the general evolution o f the interstellar 
matter (see iFalgarone fc Phillips! l)l990|) : IFalgarone et alJ 



lll995|) : ljoulain et alJ l|l998j) L For instance this is essential 
to characterize the dynamical conditions that favor the 
cooling and condensation of Hi into its molecular form. 
The same is true for the fragmentat ion/coagulation of 
dust grains l|Falgarone &; Pugetl ^995). On the theoret- 
ical side, there is a need for observational constraints for 
the ISM statistical properties. These, like the power spec- 
trum of density and velocity fluctuations in three dimen- 
sions (3D), would nicely complement the numerical simu - 
lations efforts con d ucted late l y (e.g . IPorter et alJ 1119 94') . 
IStone et al1lll998l).lMac Low! tl999)i.|Padoan &; Nordlunc 



Jl999tl. iLazarianetHuT 
lOstriker et all lj200ll) L 



200 lb . iHeitsch et alJ ll200lh 



The statistical analysis of the information contained 
in line profiles is the natural venue of such study. 
Regarding the column density fluctuations, several stud- 
ies have shown the s elf-similar properties of t he molecu- 
lar ^Falgarone et allll99ll:lBensch et all 1200 li and refer- 
ences therein) and dust emissions l)Gautier et all Il992j) 
in several regions. For Hi, the power spectrum of the 



galaxies dStanimirovic et all 


1999: Elmegreen et all 2001: 
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20011) and in the Galactic plane 


(ICrovisier & DickevL 1198.4 


Greenl Il993l iDickev et all 



Send offprint requests to: Marc-Antoine Miville-Deschenes 
Correspondence to: mamd@ias.u-psud.fr 



2001). These galactic studies were not centered on any 
particular Hi feature. Their analysis covered a depth in 
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velocity which sampled more than one spiral arm, leaving 
the necessity of studying isolated features. Furthermore, 
the statistical properties of the column density do not al- 
low to fully characterize interstellar turbulence. For in- 
compressible turbulent flows it is the velocity fluctuations 
power spectrum that follows a power law. Density fluctu- 
ations will appear self-similar in response to the velocity 
field, in certain conditions of compressibility or if the ob- 
served tracer is a passive scalar. Therefore, to study inter- 
stellar turbulence, it is important to deduce the statistical 
properties of both the velocity and density fields. 

Centroid velocity fields have been used to study the 



kinem a tics of Hn regions ll.Toncas fc RovL Il986l lO'Del 
19861 IMiville-Deschenes et all Il995l iGodbout et al 



1997) and to describe the velocity structure of star forming 
regions (see iMiesch et alJ l|l999|) and references therein) 
in which it is believed that the velocity field may be 
strongly coupled to the density field and could have a 
direct impact on t he star forming p r ocess and the ini- 
tial mass function IjScalo fc Chappelll . fl99flh . For the Hi 
gas, there has been only works done on external galax- 
ies llSpicker fc Feitzingeii Il988t IStanimirovic fc Lazarianl 
l200l|) but no observational study thriving to determine 
the 3D statistical properties of Galactic interstellar ve- 
locity fields. Recently, using fractional Bro wnian motion 
simulations, IMiville-Deschenes et alJ (j2£)03) showed that 
the power spectrum of the 3D velocity field of optically 
thin media can be determined directly from the centroid 
velocity field. 

In this paper, we propose an analysis of interferometric 
21 cm observations of the Ursa Major Galactic cirrus. The 
vast majority of the Galactic cirrus (named from their re- 
semblance to terr estrial cirrus in the IRAS all sky maps 
llLow et al.Lll984^ are composed of neutral hydrogen (Hi 
and H2) and dust; they represent the most abundant state 
(by mass) of interstellar matter in the Galaxy. We have 
selected this nearby high latitude cloud for our statistical 
analysis to minimize the importance of velocity crowd- 
ing and to study the small scale structure of the Hi. We 
have also chosen a cl oud for which the inertial ra nge of in- 
terstellar turbulence l|Tennekes fc Lumle"vlll972j) could be 
studied, meaning far from scales where energy inputs and 
dissipation occur (the vi scous dissipation scale is th ought 
to be -10 AU in the Hi - iFakarone fe Pueetl l|l995^ . The 
Urs a Major cirrus does no t show any star formation activ- 
ity l(Magnani et all Il990h and it is located at the edge of 
an expanding shell : the N orth Celestial Loop <|HeileFlll 9891 
iMeverdierks et alT.ll99lh . Therefore, it is reasonable to as- 
sume that the main mechanical energy input comes from 
the expanding shell, at scales much larger (>30 pc) than 
the field observed (—5 pc). The absence of stellar winds or 
outflows prevents the presence of energy injection at in- 
termediate scales. Finally, the observed cloud is optically 
thin at 21 cm which is a required condition to deduce the 
statistical properties of the 3D density and velocity field s 
from projected quantities (|Miville-Deschenes et al1 Eo03). 

The first three sections of the paper present the obser- 
vations (§121, the processing used to enhance the signal- 




Fig. 1. Low resolution average spectrum (from the 
Leiden/Dwingcloo survey) of the Ursa Major cirrus ob- 
served with DRAO. The spectral range covered at DRAO 
is delineated by a dotted line. 

to-noise ratio of the data (§03) and the 21 channel maps 
(§EJ)- Then we present a power spectrum analysis of the 
integrated emission map, the centroid velocity field and 
of individual channel maps (§0). This analysis allowed us 
to determine the spectral index 1 of the 3D density and 
velocity fields of the Hi gas. Our results are discussed in 

§ EI 



2. Observations 

The data were obtained using the Synthesis Radio 

Telescope of the Dominion Rad io Astrophysical 

Observatory (DRAO) ijLandecker et alll200fl. Two fields 
were observed and mosaiced. One in 1989 ij.Ioncas et all 
Il992|) and the second one in 1995. The observational 
parameters of both fields are reproduced in Table ^ The 
only difference between the two epochs was the addition 
of three new antennas to the interferometer. I n both cases 
the c alibrators were 3C147 and 3C309.1 jBaars et all 
dl977h . DRAO is also equipped with a 26m single dish 
telescope; every HI observation taken with this dish is 
calibrated using S7 (Tf, =100 K). The combination of 
the observations made by both telescopes (synthesis and 
single dish) allows the complete coverage of the u-v plane 
down to the smallest scale res olved by the synthesis 
telescope Csee lTavlor et all (|2003h where the combination 
procedure is extensively described). 

The mean 21 cm spectrum of the Ursa Major cir- 
rus, taken with a beam of 30 arcminutes with th e 
Leiden/Dwingeloo telescope ijBurton fe Hartmannl Il99^) . 
is shown in Fig^ The main Hi component is near 5 km 
s _1 but there is also an intermediate velocity component 
near -55 km s _1 , possibly related to the Perseus arm. In 
order to study the kinematics of the gas we have selected 
a channel width of 0.41 km s" 1 . The real velocity resolu- 
tion is 0.65 km s^ 1 due to filter spillover. Considering the 
number of channels of the spectrometer (see Tabled we 
have limited our observations to the main Hi component 
(see Fig^). This compromise between spectral resolution 

1 We refer to the spectral index as the exponent of the power 
spectrum assuming it can be represented by a power law. 
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Parameter 


Value 


Field 1 


ct20oo = 9ft47m 




(5 2 ooo = 70° 30' 


Field 2 


c*20oo = 9h28m 




<5 2 ooo = 70° 10' 


Channel separation 


0.412 km s -1 


Channel resolution 


0.66 km s _1 


Number of channels 


128 




-6.6 km s -1 


Vmax 


19.4 km s -1 


Angular resolution 


0.99' x 1.06' 


Field of view of each field 


2.6° 



Table 1. Log of the DRAO observations. When the two 
fields are mosaiced together, this data set contains a total 
of 107482 spectra. 



and velocity span allows us to scan more than 80% of the 
emission with a velocity resolution good enough to resolve 
the thermal broadening of the CNM. 

3. Data processing 

The preliminary steps of the reduction were done at 
DRAO (combination of the observations from the inter- 
ferometer and from the single dish, construction of the 
ch annel maps). These standard operations are described 
in lJoncas et al . (1992). The observations presented in this 
paper are characterized by a somewhat low signal-to-noise 
ratio (< 15 - see Fig. 0J. This is partly compensated in 
the statistical studies we present here by the large size of 
the sample. Once mosaiced together, the data set contains 
107482 spectra 2 . In this section we present the filtering 
method we have used to increase the signal-to-noise ratio 
of the data. 

3.1. Characterization of the noise 

The observations have been conducted such that there are 
some channels with basically no Hi emission. These chan- 
nels are used to characterize the noise. An empty channel 
of the combined Ursa Major observations (see Fig[2]- mid- 
dle) shows that the noise level is not spatially uniform. In 
the combination of the interferometric and single dish ob- 
servations we have normalized the data for the beam of 
the single dish. This operation increases the noise level 
near the edge of the field. This increase is quantified by 
computing the standard deviation of the high frequency 
fluctuations (high-pass filtering) on each spectrum of the 
data cube (see FigEl- bottom). The noise level varies from 
~ 3 to ~ 13 K from the center to the edge of the field. 
Fig |2Jtop shows the power spectrum 3 of a sub-region of 

2 As the angular resolution is ~ 1' X 1' and the pixel size is 
0.5' x 0.5', the number of independent spectra is 1/4 of that 
value (~ 26870). 

3 The technique used to compute the power spectrum is de- 
scribed in § 15.21 



001 010 1.00 




pixel 

Fig. 2. Middle: Empty channel of the combined Ursa 
Major observations. Top: Power spectrum of a sub-region 
of the empty channel (white rectangle in Middle panel). 
Bottom: Noise level at each position, estimated by the 
standard deviation of each spectrum high frequency fluc- 
tuations. 



an empty channel (see white rectangle in Fig El - middle). 
There is an increase of power at large scales (k < 0.02 
pixel -1 ) mainly due to the increase of the noise level on 
the side of the field. Part of this increase could also be 
attributed to the lower statistics of the large scale, an ef- 
fect inherent to finite size images. At small scales (k > 0.3 
pixel -1 ) the abrupt power drop is due to the beam of the 
interferometer. But in the intermediate regime the power 
spectrum is relatively flat which is an indication of the 
uniformity of the u-v plane coverage. 
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k (arcmin' 1 ) arcmin 

Fig. 3. Example of the filtering on an individual channel map. Top-right: raw data. Bottom-right: filtered data. 
Note the presence of a weak elongated structure in the south-eastern part of the field, running from (50', 80') to (90', 
40'), barely visible in the raw data and preserved by the filtering method. Left: Power spectrum of the channel map 
(inside the white rectangle area), before and after filtering. The pixel size is 30" x 30". 



3.2. Filtering 

To increase the signal-to- noise ratio of this type of spectro- 
imaging data we have developed a method based on a 
wavelet decomposition o f each channel. This meth od, an 
adaptation of the work of lStarck fc Murtaghl ljl998j) . looks 
for structures that are significant in terms of signal-to- 
noise ratio and that are also coherent in the spectral do- 
main. To optimize the filtering we have used our knowl- 
edge of the spectral and spatial variations of the noise 
level. 

First each channel cq(x, y, v) of the data cube was de- 
composed on a wavelet basis using the "a trou" algorithm: 

p 

co(x, y, v) = c p (x, y,v) + ^2 w li x -> V, v), (1) 

i=i 

where the wi(x, y, v) are the wavelet coefficients at a given 
scale I and the c v {x,y,v) are the large scale residuals. 
Then, the wavelet coefficients wi(x,y,v) are compared to 



the amplitude of the noise ni{x 1 y) at each scale and posi- 
tion. The ni(x,y) were computed as 

ni{x,y) = A n (x,y) x cr n (l) (2) 

where A n (x,y) is the noise level as a function of position 
(Fig- EJbottom) and a n (l) is the dispersion of the wavelet 
coefficients of a white noise image wit h a standard devia- 
tion of one l|Starck fc Murtaghl Il99sl) . 

If the wavelet coefficient satisfies the following crite- 
rion: 

\wi(x, y,v)\ > 1.5 x m(x, y) (3) 

it is kept, otherwise it is put to zero. Each channel is 
then reconstr ucted with the significan t wavelet coefficients 
using Eq.^ Starck fc Murtaghl <|1998|) suggest a threshold 
of 3 but we have chosen to be more conservative using a 
threshold of 1.5 to limit the removal of real weak and 
extended structures from the data. 

This filtering removes all the fluctuations that do not 
have a significant signal-to-noise ratio, but there are al- 
ways strong noise peaks that are not removed. To detect 
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Fig. 4. Typical 21 cm spectrum of the DRAO observa- 
tion, before (histogram) and after (continuous line) the 
filtering. 

these noise peaks, we use the fact that they are not corre- 
lated from one channel to the other. In practice, each spec- 
trum is checked for discontinuities larger than 8 K that go 
up and down (or down and up) in one or two channels. 
These noise peaks were removed and the missing values 
were interpolated on the spectral axis. 

An example of the filtering on a channel map is shown 
in Fig El The power spectrum of this channel map, com- 
puted before and after filtering (see Fig 01 shows that the 
noise in the raw data dominates the signal at scales smaller 
than ~ 20'. With the filtering technique applied here, we 
were able to lower the noise significantly, decreasing the 
scale range dominated by noise. In the filtered channel 
map shown here, the noise becomes important at scales 
smaller than ~ 4' (8 pixels). 

The filtering applied to the data allows the detection 
of spatially coherent structures that had a signal-to-noise 
ratio down to 1 in the raw channel maps. For instance, in 
Fig |3 we note the presence of a weak elongated structure 
in the south-eastern part of the field, running from (50', 
80') to (90', 40'), that was barely visible in the raw data. 
This very thin filament has a thickness of ~ 0.06 pc and a 
length of greater than 1 pc. The filtering has also reduced 
considerably the noise on individual spectra as shown in 
Fig. El This is of major importance for the analysis of the 
velocity field. 

4. The 21 cm data of the Ursa Major cirrus 

4.1. 21 cm channel maps 

The Ursa Major cirrus is composed of two connected parts. 
The eastern region (field 1) is predominantly atomic as it 
is detected only in the Hi line. The western region (field 



2) pr esents significant CO emission l|Pound fc Goodmanl 
Il997l). The data from field 1 has already been discussed 
bv l.Ioncas et alJ ( 19921) but the data processing techniques 
used since then have given a second life to these data. 
Figures El to show the Hi mosaic in a series of panels 
separated by 0.412 km s _1 . The velocity is indicated in 
the upper right portion of each panel. The brightness tem- 
perature in the channel map s varies from 1 K to 40 K. T he 
CO emission as observed bv lPound fc Goodmanl 1 19971) is 
present from 0.0 km s _1 to 3.71 km s _1 . It is shown in 
contours, superimposed on the Hi emission, in Figures 
tod 

The 21 cm channel maps show an intricate structure of 
interwoven elongated structures with aspect ratio up to 15. 
For the sake of simplicity, w e call them filaments in the fol - 
lowing, but, as indicated bv lLazarian k. Poeosvanl l)2000|) . 
we are aware that the structure identified in any channel 
map reflect more velocity features than density enhance- 
ments. Therefore we deliberately restrain the description 
of the structure observed in channel maps. Nevertheless it 
is interesting to point out that the western region, where 
CO emission is observed, shows generally less filamentary 
structure than the eastern part of the field where no molec- 
ular emission is detected. The detailed comparison of the 
molecular and atomic components of the cloud, which is 
clearly beyond the scope of this paper, is important to un- 
derstand the specific dynamical conditions that leads to 
the formation of long-lived molecular clouds. We leave it 
to a further study. 

4.2. Integrated emission map 

One important quantity that can be derived from 21 cm 
data cubes is the integrated emission which, for optically 
thin media, is related to the column density of Hi. The 
21 cm brightness temperature Tb(v) at a given velocity v 
and within Sv of Hi gas having a kinetic temperature 
is given by: 



T B (v)=T k {l-exp(-T(v))) 



where 



t(v) 



N(v) 



C x T k 5v ' 



(4) 



(5) 



N(v) being the gas column density at velocity v and 
within Sv, and C = 1.8 23 x 10 18 cm" 2 K" 1 km" 1 s 
llKulkarni fc HeilesL Il987h . In the optically thin regime 
(t(v) < 1), the column density can be estimated with- 
out knowing the kinetic temperature of the gas. In the 
optically thin approximation 



N(v) = C x T B {v)5v. 



(6) 



As shown in Fig. |S1 the effect of opacity on the de- 
termination of the column density is stronger for cold 
gas and for high brightness temperatures. For warm Hi 
(Tfc > 500K), the opacity is negligible and the column 
density is well approximated by Eq.El On the other hand, 
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Fig. 5. Ursa Major. Channel maps of the 21 cm line emission be tween -1.65 and 2.88 km s 1 . The contours (1.5, 2 
and 3 K) represent the CO emission of lPound & Goodmanl (|l997|) . 



for typical CNM kinetic temperature (~100 K) and for 
brightness temperature of ~40 K (the maximum value in 
our channel maps) , the column density in a given channel 
map is underestimated by ~20%. But the 21 cm profiles 
are relatively broad and most of the brightness tempera- 
ture values are below 20 K. Therefore, if the kinetic tem- 
perature of the Hi in Ursa Major is ^100 K, we think 
that the column density is probably underestimated by 
only ~ 10% in the brightest region of the cloud, which 
confirms that we are in the optically thin regime. 

Instrumental noise also affects the determination of the 
column density. Since the uncertainty on the column den- 
sity increases with the square root of the number of inde- 
pendent channels used in the summation, it is common to 
define a spectral window within which the column density 



is computed. This way noise contribution from the empty 
channels is avoided. However for these observations the 
baseline is very short (see Fig. and varies in length at 
both ends as a function of the Dopplcr shifts of the line 
emission. Therefore all channels were summed. 



The integrated emission map of the Ursa Major field 
is shown in Fig. [5] This map is characterized by a set 
of filaments, mainly oriented west-east. The integrated 
emission peaks at W « 330 K km s _1 (corresponding 
to Nh ~ 6 x 10 20 cm~ 2 ) in the eastern part of the main 
filament. The mean integrated emission value in the field 



is W « 190 K km s" 1 (N H ~ 3.45 x 10 20 cm 
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Fig. 6. Ursa Major. Channel maps of the 21 cm line emission b etween 3.30 and 7.83 km s 1 . The contours (1.5, 2 
and 3 K) represent the CO emission of lPound fc Goodman! (|l997|) . 



4.3. Centroid velocity field 

The centroid velocity field C(a,5) of a position-position- 
velocity (PPV) data cube Te(a,(5, v) is given by the fol- 
lowing expression: 



(7) 



nl c\ J2v vx T B (a,8,v) _ 1 
C a, <5) = — ^— km s . 

For optically thin clouds, the centroid velocity field is ex- 
actly the mean ra dial v e locity of the gas on the line of 
sight. As shown bv lPetvl 1^999), noise significantly affects 
the determination of velocity centroids. In fact, the er- 
ror on the velocity centroid determination depends on the 
signal-to-noise ratio of the spectrum and on the number 
of channels used. To minimize the effect of noise, one may 
compute the centroid velocity on each spectrum by se- 



lecting a spectral window on each spectrum where all the 
brightness temperatures are greater than twice the noise 
level. But practically, the wavelet filtering we have applied 
significantly reduces the noise level of the data leaving very 
few noise peaks that could contaminate the centroid veloc- 
ity estimate. Therefore, and as there are very few channels 
with no emission, all the channels were used to compute 
the centroid velocity of each spectra. The centroid velocity 
field of the Ursa Major cirrus is shown in Fig. 

The centroid velocity field of the Ursa Major cirrus 
is characterized by a set of elongated structures. In fact, 
the two main elongated structures seen in the integrated 
emission map and in channel maps also stand out in the 
centroid velocity field. Velocity fluctuations are observed 
at all scales, up to the largest one where there is an east- 
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Fig. 7. Ursa Major. Channel maps of the 21 cm line emission between 8.24 and 12.77 km s 



west velocity gradient of ~ 1.5 km s _1 pc _1 at a 2° scale 
(considering a distance of 100 pc for the cirrus). But it is 
at small scales that the velocity gradients are the largest. 
Across the brightest filament of the cirrus (near position 
(50', 80') in Fig. HUI) we observe transverse velocity gradi- 
ents betw een ~ 5 and 10 km s" 1 PC" 1 , a t a 20' scale (also 
shown bv lMiville-Deschenes etafl feOOtft ). 

5. Power spectrum analysis 

5.1. Determination of the density and velocity power 
spectra in three dimensions 

One important aspect of the study of interstellar 
turbulence is to determine the statistical properties 
of the density and velocity fields in three dimen- 
sions (3D). To do so one must understand the ef- 



fect of opacity and of projection on the observed 
quanti ties. Regarding density fluctuations, many au- 
thors jS tutzki e t all Il99 8t lLazarian k. Poeosvanl [2000; 
iGoldmanl l2000t iMiville-Deschenes et all l2003|) showed 
that the power spectrum of the 3D density field can be 
determined directly from the power spectrum of the inte- 
grated emission map under two conditions: 1) the observed 
medium must be optically thin and 2) the spatial scales 
observed on the plane of the sky must be smaller or equal 
to the spatial depth of the line of sight. Under these condi- 
tions, the power spectrum of the integrated emission map 
is proportional 4 to the power spectrum of the 3D density 
field. 



4 It has the same shape - 
normalisation. 



or spectral index - but not the same 
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Fig. 8. Error done on the column density estimate using 
the optically thin approximation (Eq. EJ) for Hi gas with a 
kinetic temperature of 50, 100 and 500 K. Bottom: The 
column density (per velocity interval) determined using 
the optically thin approximation as a function of bright- 
ness temperature is given by the dashed line. The three 
continuous lines indicate the real column density (per ve- 
locity interval) of the gas as a function of brightness tem- 
perature. Top: Fraction of the real column density that is 
missed in the optically thin approximation as a function 
of the brightness temperature at a given velocity. 



The determination of the 3D velocity field power spec- 
trum is less obvious. lLazarian &: Poeosvanl l)2000|) pro- 
posed a method called Velocity Channel Analysis (VCA) 
that uses the variation of the spectral index of veloc- 
ity slices of increasing thickne ss. But as pointed out by 
IMiville-Deschenes et al.l (|^003) and as we will show here, 
this method is ver y difficult to apply on rea l obser vations. 
On the other hand. IMiville-Deschenes et alJ l)2003[) showed 
that centroid velocity fields can be used to determine the 
power spectrum of the 3D velocity field. Using Fractional 
Brownian Motion simulations 5 these authors studied the 
impact of density and velocity fluctuations on the power 
spectrum of centroid velocity fields. They showed that 



5 Fractional Brownian motion simulations produce self- 
similar Gaussian random images. 



density fluctuations have a very limited effect and that the 
power spectrum of the centroid velocity field is the same 
as for the 3D velocity field, assuming that the medium 
observed is optically thin and that the scales observed on 
the sky are smaller than the depth of the line of sight. 

It is important to estimate if these results can be ap- 
plied to the specific physical conditions of the Ursa Major 
cirrus. First, we have shown that opacity effects are limited 
for the brightness temperature observed here (see Fig. |SJ . 
Second, as the cloud is at 100 pc, the largest transverse 
scale in the field is only ~7 pc, which is very likely to 
be smaller than the depth of the Hi froth in that region 
(probably several tens of pc if we look at large scale obser- 
vations - see Fig. 1120 . Therefore, in the following we an- 
alyze the power spectrum of the integrated emission map 
and of the centroid velocity field of the Ursa Major cirrus 
in the framework of IMiville-Deschenes et al" 1 ll2003h . 

5.2. Integrated emission map and centroid velocity 
field 

For the Ursa Major observations, we have limited the 
power spectrum computation to a rectangular region 
(354 x 198 pixels) where all the points have been observed. 
This sub-region is delimited by a white rectangle in fig- 
ures and ^]] The power spectrum is the sum of the 
square of the real and imaginary parts of the two dimen- 
sional FFT, azimuthally averaged over each wavenumber 
k = \lk% + ky. As stated bv lBensch et al.ll|200l|) . edge ef- 
fects may introduce systematic bias in the determination 
of the power spectrum index using Fast Fourier Transform 
algorithms. When the power spectrum is described by a 
power law P(k) oc k@, the error on the power spectrum in- 
dex depends on (3 itself; it becomes important for [3 < —3. 
To reduce this effect, we have apodised the edge of the 
image with a cosine function. The size of the apodised 
region is equal to 3% of the linear size of the image. 
We have tested this met hod on non-periodic Fractiona l 
Brownian Motion images l|Miville-Deschenes et all 2003) 
and we have found that the p value is recovered with a 
typical error of 0.1 for —4 < (3 < —2. 

The power spectra (raw and filtered data) of the inte- 
grated emission map and centroid velocity field are shown 
in Fig. [§] and ^3 For both the integrated emission map 
and the centroid velocity field the power spectrum of the 
filtered data is well described by a power law oc fc- 3 - 6 ±°- 2 
for wavenumbers 0.01 < k < 0.2 arcmin -1 . We have to re- 
strict the slope determination to a relatively narrow scale 
range for two reasons. First, at large scale, the ID power 
spectrum points have a large uncertainty as they are the 
average of a small number of samples in the 2D FFT. 
Second, at small scales, the power spectrum is affected by 
the beam, the noise and the non-periodicity of the image 
mentioned earlier. The uncertainty on the spectral index 
we give in this paper (±0.2) reflects the dispersion of the 
Pk values around the power law but mainly the uncer- 
tainty related to the apodization, beam and noise (see 
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Fig. 9. Hi Integrated emission of the Ursa Major cirrus Right-top: raw data. Right-bottom: filtered data. Left: 
Power spectrum of the integrated emission map of the raw and filtered data. The dotted line represents a power law 
with a spectral index of -3.6. 



above). This uncertainty value has been validated using 
the fractional Brownian motion simulations and realistic 
noise and beam shape. 

As opacity effects are negligeable at 21 cm and the 
depth of the cirrus is most probably larger than the largest 
transverse scale, we estimate that the 3D density and ve- 
locity fields of the Ursa Major cirrus are characterized by 
a spectral index j v = j n = —3.6 ± 0.2. Here we consider 
that the energy driving the turbulent motions is injected 
at much larger scales than the size of the observed fields 
and that the behavior of the velocity field at all scales 
is inherent to the turbulent cascade. It is interesting to 
point out that the velocity spectral index is compatible 
with the value expected (-11/3) for Kolmogorov turbu- 
lence l|Kolmogorovt Il94l|) . 

5.3. Velocity channel analysis 

As discussed earlier, lLazarian fc Pogosvanl l|2000j) pro- 
posed a method to determine the spectral index of the 
3D density and velocity fields using PPV data cubes 
of optically thin media. This method has been used 



by a few authors to determine the 3D velocity spec- 
tral i ndex j v of Hi in the Ga lactic plane i|Dic key e\ all 
| 200lh or in other galaxies |s tanimiroyic el; al l llflflflt 
IStanimirovic fc LazariaiJ, l200ll lElmeereen et al.j l200l|l . 
The VCA method is based on t he analysis of veloc- 
ity sli ces of increasing width AV. lLazarian fc Pogosvanl 
(200(3) claim that the spectral index of velocity slices 
should reach two asymptotic states at small and large AV, 
from which j v could be determined. 

We have tested the VCA method on our observations. 
Fig. El shows the variation of the spectral index of veloc- 
ity slices of increasing width AV. Here a velocity slice of 
width AV is the sum of the N = AV/ Su brightest con- 
secutive channel maps, where 5u — 0.412 km s . The 
spectral index of each velocity slice was computed; it is 
the slope of the power spectrum between 0.01 < k < 0.1 
arcmin -1 . As expected, at large AV the curve in Fig. 1111 
converges to the spectral index of the integrated emis- 
sion (^-3.6). At the other end, the curve reaches an- 
other asymptotic state to the value (~-3.1) found for the 
brightnest channel map (see Fig. . 
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Fig. 10. Hi centroid velocity field of the Ursa Major cirrus Right-top: raw data. Right-bottom: filtered data. Left: 
Power spectrum of the centroid velocity field of the raw and filtered data. The dotted line represents a power law with 
a spectral index of -3.6. 



But here it is important to estimate the effect of the 
gas temperature on this curve. The vertical dashed line in 
Fig.llllshows the effective spectral resolution of our obser- 
vations Sv e ff w y/ 'Su 2 + 2.16fcsT/m = 1.7 km s _1 (where 
m is the mean atomic mass and if we consider typical C NM 
gas at T = 150 K - see lMiville-Deschenes et al.ll|2003|) for 
details). The plateau seen at small AV in Fig. 1111 ia sim- 
ply the result of the gas temperature that washed out 
all velocity fluctuations in channel maps. Therefore, the 
asymptotic state at small A V is not the one expected by 
lLazarian fc PogosvarJ (2000) and cannot be used to deter- 
mine Like lMiville-Deschenes et alJ (|2003j) we conclude 
that while it could be used for large scale observations (i.e. 
external galaxies), the VCA technique is not suited to de- 
termine the 3D velocity spectral index for local HI where 
the velocity dispersion is commensurate with the thermal 
broadening on a broad range of scales. 

5.4. Combination with Leiden observations 

To obtain a more global view of the Hi power spec- 
tra (integrated emission and centroid velocity) we com- 



bined the power spectra computed on the Ursa Major 
DRAO observations with the power spectra computed on 
the Leiden/Dwingeloo observations of the whole North 
Celestial Loop region (which includes the Ursa Major cir- 
rus) as observed by the Leiden/Dwingeloo survey (see 
Figures IT21 and ILUl. The angular scales of each observa- 
tion complement very well to give the power spectrum 
of the Hi integrated emission and of the centroid veloc- 
ity fields on almost three orders of magnitudes in scales, 
from 2 arcminutes to 16 degrees (see Figures ITU and ITS')) . 
One outstanding result of this analysis is the fact that 
both power spectra are compatible with a single slope of 
— 3.6 ± 0.2 over the whole range of scales. The decrease 
of power at small scales (k > 0.4 arcmin^ 1 ) is due to the 
domination of the noise at these scales. All the power at 
these scales was removed by the wavelet filtering. At the 
other end of the power spectrum, the flattening may not 
be significant as the statistics at large scale are low. If we 
consider that most of the gas in the North Cele stial Loop 
is loc al, at an approximate distance of 100 pc l)Penprasel 
Il993h . and the opacity effects are limited in the region, 
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Fig. 11. Power spectrum index of velocity slices as a func- 
tion of the slice width AV. Velocity slices are built by 
adding velocity channel maps from the the brightest to 
the faintest. 

the local Hi is self-similar on spatial scales from 0.5 to 25 
pc. 

5.5. Density-velocity correlation 

The structure of the centroid velocity field of the DRAO 
fields shows many similarities with the column density 
map. In both maps, the main filaments of the cirrus are 
visible, indicating that the density and velocity of the Hi 
show some correlation. This is also true at larger scales as 
one can see in Figures IT21 and IT51 where we show the inte- 
grated emission and the centroid velocity field of the North 
Celestial Loop (Leiden/Dwingeloo data). At the scale of 
the loop, the main density structures seen in the integrated 
emission map are correlated with structure in the centroid 
velocity field. 

To estimate the correlation level between the column 
density and the velocity field, we have computed the cross- 
correlation defined as: 



Ccr 



(N H tot V c ) - (AW) (V c ) 



Cat cry 



(8) 



where erjv and ay are the standard deviation of the col- 
umn density and centroid velocity fields respectively. The 
cross-correlation is 0.44 for the Ursa Major field and 0.25 
for the North Celestial Loop region, indicating a moder- 
ate correlation between the centroid velocity and column 
density field. 



6. Discussion 

As shown by 



lLazarian fc Pogosvanl fcOOOTl and 
IMiville-Deschenes et alJ l|2003|) . the interpretation of 
the power spectrum of individual channel maps is not 
unique since brightness fluctuations in a channel map 
are the combination of velocity, density and temperature 



fluctuations. On the other hand, for optically thin and 
isolated regions, and at scales that are smaller than the 
depth of the observed region, the power spectrum of the 
Hi integrated emission is exactly the power spectrum of 
the 3D density field. Furthermor e, under these conditions 
and in the case of Gaussian fields. IMiville-Deschenes et alJ 
(2003) show that the power spectrum of the centroid 
velocity field is exactly the power spectrum of the 3D 
velocity field. 

With these ideas in mind, it is interesting to 
compare our results (j n = —3.6 ± 0.2 and j v = 
—3.6 ± 0.2) with previous power spectrum analysis of 
21 cm observations, and try to draw a coherent pic- 
ture of the Hi density and ve l ocity st atistica l prot 



erties. First, 



Croyisier fc Dickevl (119831): lOreenl ijlQc 
iDeshnande et all (|2000lk iDickev et all l(200ll) all found a 
spectral index of ~ —2.9 ±0.2 for individual channel maps 
of Galactic plane regions with significant emission, in ac- 
cordance with our values for individual channels in the line 
center (see Figures l3l and ITT|l . But, as it has been said be- 
fore, the contribution of the velocity, density and temper- 
ature fluctuations to the structure in individual channel 
maps makes these spectral indexes ver y difficult to inter- 
pret. Aware of this. lDickev et al.l l|200l|) also looked at the 
spectral index of the integrated emission which, in their 
most out-of-plane region (b = 1.3°), is ~ —4. However, 
one must keep in mind that the power spectrum of an 
integrated emission map of Galactic plane regions could 
be affected by non-turbulent large scale structure and 
by significant Hi self-ab sorption. Recent hi gh-resolution 
(1') 21 cm observations I Gibson et all EoOff) have shown 
the presence of several Hi self-absorption regions in the 
Galactic plane, that could not be detected in previous low- 
resolution 21 cm surveys. As self-absorption seems to be 
more important at small-scale, it could reduce significantly 
the amplitude of the small-scale Hi emission and produce 
an artificial steepeni ng of the power spectru m. Lets also 
mention the work of iDeshpande et all l|2000|) who looked 
at the power spectrum of the velocity-integrated Hi opti- 
cal depth, using absorption measurements in the direction 
of Cassiopeia A. They found a relatively high spectral in- 
dex (—2.75 ± 0.25) but, as they pointed out, the optical 
depth depends both on the column density and the gas 
temperature which makes the relation with the 3D den- 
sity power spectrum not straightforward. 

Recently there have also been some power spectrum 
analys is of 21 cm emission in externa l galaxies. For in- 
stance Istanimirovic fc Lazarianl l|200ll) found j n = —3.3 
for the Hi integrated emission in the Small Magellanic 
Cloud (scales from 30 pc to4kpc) . In the Large Magellanic 
Cloud. lElmegreen et alJ ( 200ll) found that the power spec- 
trum of the Hi integrated emission is well described by two 
power laws of indexes -3.66 and -2.66, on spatial scales be- 
tween 20 pc and 9 kpc. The spectral index change, which 
occurs near L ~ 100 pc, can be explained by the fact that 
the scales on the plane of the sky are larger than the depth 
of the medium for L > 100 pc. In this case, the steepening 
of the slope observed at small scales reflects a change in 
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Fig. 12. Integrated emission of the 21 cm 
Leiden/Dwingeloo survey in the direction of the North 
Celestial Loop. The Ursa Major fields observed with the 
DRAO interferometer are indicated by boxes. 



the topology of the medium investigated from a 2D (large 
sca les) to a 3D (small scale s ) syst em, an effect described 
by iMiville-Deschenes et al.l ll2003h. The de nsity spectral 
index measured by lElmegreen et al in the LMC 

is in accordance with our j n value, even if their smallest 
scale corresponds approxi mately to our largest sca le. On 
the other hand, contrary to lElmegreen et alJfcOOlh we do 
not observe a break in the slope of the power spectrum 
which indicates that the scale range we trace (from 0.1 to 
25 pc) is smaller than the typical depth of the medium. 
This is in agreement wit h the fact that the sca l e heig ht of 
Galactic Hi is ~ 100 pc l|Lockman &: Gehmanl . fl99lj) . 

The diversity of Hi spectral indexes found in the litter- 
ature reveals the importance of gas temperature, depth, 
energy injection and opacity effects on the power spec- 
trum of 21 cm emission and how difficult it is to determine 
accurately the 3D statistical properties of interstellar tur- 
bulence. Unlike most of the Hi studies done so far, 21 cm 
observations of high latitude regions like the one presented 
in this work offer the most favorable physical conditions to 
determine the 3D density and velocity power spectrum of 
interstellar turbulence. Such regions are almost unaffected 
by self-absorption and their density and velocity fields are 
dominated by turbulent motions, which is not always the 
case for Galactic plane regions and external galaxies where 
energy injection at several scales by supernovae, stellar 
winds or outflows might affect the power spectrum. The 
work we present in this paper is the first study of the in- 
ertial range of interstellar turbulence, far from energy in- 
jection and dissipation scales, and where opacity and pro- 
jections effects are controled. It is the first time that the 
3D density and velocity power spectra, essential quantities 
for the understanding of interstellar turbulence, are deter- 



mined for the Galactic interstellar medium. Furthermore, 
the dynamical range of scales over which they are mea- 
sured (more than two orders of magnitude) is larger than 
for all previous studies. The existence of these two scaling 
laws is therefore a robust result. 




a m (hr> 



Fig. 13. Centroid velocity field of North Celestial Loop 
(Leiden/Dwingeloo data). See Fig. 1121 

It is important to point out that our results (7„ = 
— 3.6±0.2 and j v — — 3.6±0.2) agree with the Kolmogorov 
prediction for incompressible turbulence (-11/3). It may 
indicate that compressibility and magnetic field do not 
affect significantly the energy cascade of interstellar tur- 
bu lence. Our results are in a ccordance with the work 
of IGoldreich fe Sridharl l|l995l) (see also the review of 
ICho et all ((2003)) who predicts a Kolmogorov type en- 
ergy spectrum for an incompressible turb ulent magne- 
tized fluid. But IGoldreich fc Sridharl Jl995) also predicts 
an anisotropy of the velocity and density fluctuations that 
increases at small scales. It is interesting to point out here 
that the density and velocity structures of the Ursa Major 
cirrus are clearly anisotropic; filaments with aspect ratio of 
15 are observed. One important step in characterising in- 
terstellar turbulence would be to quantify this anisotropy 
and its variation with scale. 

It is also interesting to note that the spectral indexes 
of the Hi are significantly different from what is observed 
with molecular or dust emission. Several studies of inte- 
grated CO emission shows spectral inde xes around -2.5 
and -2.8. For instance |r3ensch et alJ (|200ll ) found relatively 
uniform spectral index values (—2.8 < (3 < —2.5) for a set 
of five molecular clouds observed with tracers of different 
optical depth ( 12 CO and 13 CO J = 1 — » 0). The work of 
iGautier et alJ l)l992j) on 100 /im dust emission also shows 
a relatively shallow slope j3 ~ —3 for non-starforming 
regions. These results seem to indicate more small scale 
structures than what is observed in Hi. However molecular 
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Fig. 14. Power spectrum of the integrated emission of 
the Ursa Major / North Celestial Loop region. The dot- 
ted line defines the scale regimes of the DRAO and 
Lcidcn/Dwingeloo data (see Fig. I12j) . The dashed line rep- 
resents a slope of -3.6. 
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Fig. 15. Power spectrum of the centroid velocity field 
of the Ursa Major / North Celestial Loop region. The 
dotted line defines the scale regimes of the DRAO and 
Lcidcn/Dwingeloo data (see Fig. ll3fl . The dashed line rep- 
resents a slope of -3.6. 



emission is, contrary to Hi, strongly affected by opacity 
and heating variations within the cloud which may sig- 
nificantly modify the power spectrum. Similarly the dust 
thermal e mission may vary locally accor ding to the grain 
properties ([Miville-Deschenes et all 120021) . to the extinc- 
tion of the stellar radiation or the abundance of H2. This 
apparent difference between the atomic and molecular gas 
structure clearly deserves further analysis. 



7. Conclusion 

We have presented a power spectrum analysis of 21 cm 
interferometric observations of the Ursa Major / North 
Celestial Loop region. These observations of an isolated 
high-latitude region, with no star-formation activity, al- 
low to study the statistical properties of the inertial range 
cascade of interstellar turbulence in the Hi. From the anal- 
ysis of the integrated emission, of velocity channels and of 
the centroid velocity field of the Hi emission with a few 
assumptions on the depth of the emitting medium and its 
transparency, we were able to estimate the spectral in- 
dexes of the density and velocity fields in 3-D. For both 
fields we find a spectral index of — 3.6±0.2. Last the veloc- 
ity field and the integrated emission map show moderate 
correlation with a value of at most 0.4 in the Ursa Major 
field. 
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